Elliptic flow at varying energy heavy ion collisions: Partonic versus hadronic dynamics 
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We examine whether the breakdown in elliptic flow quark number scaling observed at the Rel- 
ativistic Heavy Ion Collider (RHIC) energy scan is related to the turning off of deconfinement by 
testing the hypothesis that hydrodynamics and parton coalescence always apply, but are obscured, 
at lower energies, by variations in the widths of quark and anti-quark rapidity distribution. We 
find that this effect is enough to spoil quark number scaling in elliptic flow. A lack of scaling in 
data therefore does not signal the absence of partonic degrees of freedom and hadronization by 
coalescence. In a coalescing partonic fluid, however, elliptic flow of anti-baryons should be greater 
than that of baryons, since anti-baryons contain a greater admixture of partons from the highly 
flowing mid-rapidity region. Intriguingly, purely hadronic dynamics has a similar dependence of 
baryon-anti-baryons elliptic flow as purely partonic dynamics, again because anti-baryons tend to 
come from regions where the deviation of the system from hydrodynamic behavior is at its smallest. 
The opposite trend observed in experiment is therefore an indication that we might be misunder- 
standing the origin of elliptic flow. We finish by discussing possible explanations of this, and suggest 
experimental measurements capable of clarifying the situation. 
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PACS numbers: 13.87.-a, 12. 38. Aw, 25.75.-q, 24.60.-k 



I. INTRODUCTION 

The discovery, at the relativistic heavy ion collider 
(RHIC), of a "perfect fluid" hadronizing via "quark co- 
alescence" [l|-l7| has aroused a lot of experimental and 
theoretical attention. The main experimental evidence 
to support perfect fluidity is the observation that ellip- 
tic flow V2, defined as the 2nd Fourier coefficient in the 
particle spectrum wrt the reaction plane <f> rp 
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(1 + 2v n cos (n(0 - 4>rp))) (1) 



approaches the value expected from hydrodynamics [fj. 
The further scaling of mesonic and baryonic V2 is sugges- 
tive of what is seen by the most naive quark coalescence: 
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Provided the Wigner functions [8|, Q W(...) are 5- 
functions in position and momentum, i>2 ar!/on (p/3)/3 = 
v 7 2 neson {p/2)/2 [n| , a relation that holds surprisingly 
well at the top RHIC energy of 200 GcV per nucleon 
[Tj-Q , even if signs of finite width Wigner functions have 
been identified 0, [ll| . The most common interpretation 
of this is that the deg rees of freedom carrying elliptic 
flow are partonic |7H13j . Coalescence of partonic degrees 



of freedom is also hinted by the high pt baryon distri- 
bution [1, El i an d can account for 2-particle correlations 
structures such as the near-side "ridge" [l2[ and the far- 
side "cone" [13]. It has been found to significantly affect 
also the nuclear modification factor and the elliptic flow 
of heavy quarks [H, EH . 

It is, however, far from clear how the flow data and 
its "partonic ideal fluid" interpretation is related to the 
onset of deconfinement. In the confined phase, the cou- 
pling constant between mesons goes as A/iV 2 [17| while 
the density of degrees of freedom, impacting both vis- 
cosity r\ and entropy density s, goes as N®. Hence, 
r]/s ~ N% . In the deconfined phase the density of degrees 
of freedom goes as AT?, hence rj/s ~ N® independently 
of the coupling constant. t]/s should respectively go to 
a constant for strongly interacting theories [lH or go as 
~ (A 2 In A) -1 in the perturbative regime fl9ll. both ~ N®. 
Hence, Simple number of colors scaling [131 shows that 
at deconfinement rj/s should jump by an order of mag- 
nitude. If elliptic flow is indeed hydrodynamically gen- 
erated, a corresponding jump should be observed when 
the initial temperature is about the critical temperature 
for deconfinement [2(| [2l|. No such jumps in the ellip- 
tic flow observable are apparent in experimental data, 
whether varied in center of mass energy, centrality or ra- 
pidity SHHH. 

Seeing how the v-i observable "turns off" is indeed one 
of the objectives of making lower energy measurements 
with the latest detectors [23l - l26| . Recently, the RHIC 
experimental energy scan started observing a system- 
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FIG. 1. (color online) dN/dy of quarks and anti-quarks. For 
anti-quarks we plot the 3 widths employed in our studies, 
a y = 0.25,0.5, 1. 



atic breakdown of coalescence in yfs (not seen, so far, in 
rapidity [2?]]) which has been interpreted as signifying 
the "turning off" of partonic degrees of freedom [28|, [29| . 
However, as pointed out in [3(|, such conclusion is prema- 
ture: even if quark coalescence persists, high baryochem- 
ical potential and small rapidity intervals in the collision 
region might be enough to make it not apparent when vi 
of mesons and baryons is considered. 

In this work we explore the consequences of the rapid- 
ity dependence of both dN/dy and V2{y) on the quark 
number scaling and on the baryon and anti-baryons split- 
ting of the elliptic flow, pointing out that recent observa- 
tions of the RHIC experiments at low energy are not ex- 
plicable neither by a pure hydrodynamical+coalcscence 
model nor by a pure hadronic model. 



II. MODEL 

Our model is simplified to capture the main issues we 
want to deal with. We assume all v 2 is produced in a 
closc-to-idcal hydrodynamic stage at all energies. 

The main consequence of an ideal hydrodynamics stage 
for our purposes is that, locally in rapidity and co- moving 
time, flow is species independent. For configuration 
space coordinates x, while abundances of flavors can vary, 



««(£) = u«(f) = «'(f) 



(3) 



This is simply a consequence of the fact that in ideal 
fluids any conserved current J M cx u^, and deviations 
from this relation are invariably dissipative (3lj . This 
is independent of the initial condition for hydrodynamic 
evolution allowing to draw qualitative conclusions from 
an admittedly simplified schematic model. 

Hence, before hadronization, provided flavor content is 
homogeneous in the transverse plane, for quarks emitted 



at rapidity y and transverse momentum px we have, up 
to mass corrections 



,,i' s 



(pr,y) = vf s (p T ,y), (pT) q ' s (y) = (pt) 11 ' 3 (y) (4) 



The other important ingredient to define the partonic 
stage is the y and pt dependence of quark and antiquark 
distribution function. We assume a Gaussian antiquark 
distribution (Fig. [T]) 
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f(PT, 0) 



(5) 



where a 2 , is approximately ~ Ins, growing from ~ 
]n(s/2m pr oton) to ~ ln(s/m pro t on ) between the Landau 
and Bjorken limits (32l. l33j. 

For the quark distribution, we assume a flat distribu- 
tion for a rapidity region \y\ < 1 meaning the probabil- 
ity of finding a valence quark is approximately invariant 
with rapidity. This appears quite reasonable at least if 
the beam energy is y/s > 5GcV, 
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= N°f( PT ,< 



\d 2 p T dy / 

The normalization can be fixed by the p/p ratio 



(dN°/dy)/(dN° q /dy) ~ [(dN ¥ /dy)/(dN p /dy)} 



1/3 



(6) 



(7) 



experimentally decreasing with yfs, as more entropy is 
carried by pre-existing valence quarks rather than created 
qq pairs. Strange quarks and antiquarks distribution is 
the same as antiquarks, since they are produced from 
zero in the hot medium. 

While these are undoubtedly simplified assumptions, 
more complicated scenarios [34[ would, provided they 
model the experimentally observed limiting fragmenta- 
tion of t>2 and dN/dy [2l[ and Eq. [3J give qualitatively 
similar results to those we present. This is because in any 
such system the bulk of antiquarks and strange quarks 
would come from a central plateau, whose width Ay < 1 
grows very slowly with yfs [H, [2l[ . The central plateau 
is also the region closest to ideal hydrodynamics, at all 
energies. 

In accord with the hydrodynamic hypothesis for the 
quark phase, the momentum distribution is independent 
of the quark flavor (q,q, s,s have the same {pt) and t^)- 
Hence, a similar f(pr, 4 1 ) describes both q and q 

f(p T , 0) = e-^ [1 + 2v p 2 (p T , y) cos (20)] (8) 

and T = 170MeV. To fit experimental data [HlHIl, the 
quark (antiquark and strange quark) v 2 distribution for 
partons is also Gaussian: 



vlipT,y) = v 2 (y = 0,p T )e "y 



(9) 



v^iy = 0,Pt) is an empirical function compatible with 
experimental data. Its integral is fixed by the experimen- 
tal constraint [U v 2 (y = 0)/e - (1/S) (dN/dy), where 
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FIG. 2. (color online) V2/n q vs pr/riq of quarks, mesons 
baryons and anti-baryons in the \y\ < 0.5 rapidity window and 
for two different cr values: Top panel (a) assumes a y — 0.25 
and bottom panel (b) has a y = 0.5. 



S is the overlap area of the collision S ~ (1 — e)irR 2 . 
V2(y > 0) is assumed, as seems to be the case in experi- 
mental data [l|, [U, , to have the same a y as the "hot" 
medium (antiquark and strange quark distributions). 

Note that even if V2(y,PT) is the same for all partons 
locally, this is not true for V2{pt) integrated over a finite 
y region. For example, in Fig. 2, we show V2(pt)||j,|<o.5 
for quarks (up triangles) and antiquarks/strange quarks 
(down triangles) for two different widths. 

This difference is due to the different dN/dy distri- 
butions for quarks and antiquarks together with the ra- 
pidity dependence of V2- Anti-quarks have the largest 
dN/dy from y ~ region because at finite y dN/dy has 
a Gaussian-like tail 

-y 2 

. Jly.dye ^v 2 ( PT ,y) 

V2,q{PT)\\y\<y> = ~ 2 (10) 



Since v 2 is maximum at y ~ 0, this generally leads to 
V2,q(pr) < V2,q(pr)- Since the Wigner function is also 
approximately a Gaussian, a similar mechanism acts also 
in the coalescence process leading to a smaller flow for 



FIG. 3. (color online) V2 of baryons and antibaryons calcu- 
lated in the rapidity region \y\ < 0.5. The antiquark and 
strange quark distribution is a Gaussian with width in rapid- 
ity of, respectively, a y — 0.25 (top panel (a) ) and a y = 0.5 
(bottom panel (b) ). 



mesons with both an q and an s quark, such as K~ , with 
respect to those having a qq or qs. 

To study the v 2 of hadrons in the px range [0 — 5] GeV 
we assume coalescence-type hadronization that for the 
case of mesons is given by 

n m =c m j n (p da ^ s (Pi ~ m ^ 

i—a.b 

fi{ri,Pi)W (r| ) • (ii) 

The relative phase space coordinates r = ri } — r a and 
q = pb — p a are the four- vector relative space-time dis- 
tance and energy-momentum, da is a volume element 
of a space-like hyper-surface. The hyper-surface of co- 
alescing partons is fixed through the condition of equal 
longitudinal proper time r = Vi 2 — z 2 [Hj]. In the 
coalescence integral in Eq. (|TTj) we consider the full phase 
space overlap of the coalescing particles, with the advan- 
tage of avoiding some of the more restrictive approxima- 
tions employed by other coalescence models. In particu- 
lar the extended 3D Wigner function avoids any collinear 
approximations in momentum. 



The hadron Wigner function for light quarks used is 
a simple product of spheres in position and momentum 
space 
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(12) 



The radii A r and A p in the Wigner formalism obey the 
relation A p = A" 1 , motivated by the uncertainty princi- 
ple. Similar expression can be derived for 3-quarks coa- 
lescence. The parameter A p is taken to be different for 
baryons and mesons and is of the order of the Fermi mo- 
mentum, A^ = 0.23 GeV and A^ = 0.35 GeV . This has 
been shown to account for both the it, K, p,p, A spectra 
at RHIC energy and quark number scaling elliptic flow 
including its small baryon/meson breaking [Til l35l . |36| . 
The condition for the 4- vector relative momentum can 
be written as: 



(PTm - PT,b) 2 - 2mx, a m T ,b cosh(y a - y b ) 



T.a 



l T.b 



(13) 



in the non-relativistic limit \p\ << mi this reduces to 
the condition (p a ~ p b ) 2 < A p , i.e. that the relative mo- 
mentum, and hence also the difference in rapidity, of the 
coalescing quarks is limited to the width of the hadron 
wave function. 

In the ultra-relativistic limit \p\ >> rrii Eq. (fT3"|) re- 
duces to: 

2p Tta pT,b[ cos ( ( f ) a - <j>b) + cosh(y a - y b )] < A 2 p (14) 

which implies that at large px , due to the Lorentz boost 
from the fireball frame to the coalescing quark frame, 
particles with relatively large Ap = \px, a — Pr,b\ can coa- 
lesce. This effect [ll[ follows from relativistic kinematics 
and leads in our case to an increasing difference between 
particle and anti-particle V2(j>t) as a function of px- 

While in previous work the elliptic flow of coalescence 
was studied considering flat distribution in rapidity, we 
show that in presence of a finite Wigner function in Eq. 
[2] and a y— dependent anti-quark density the naive quark 
number scaling breaks down and a difference in baryon- 
antibaryon vi appears even if ideal partonic hydrodynam- 
ics and coalescence occur at all energies. While such a 
scenario is usually assumed not to apply for lower en- 
ergies (this is very different from "demonstratively fal- 
sified"), it is worthwhile to estimate the magnitude of 
these effects. 



III. RESULTS 

The meson and baryon i^'s are summarized in Fig. [5J 
As can be seen, the v n {px /n q ) jn q naive scaling is con- 
siderably broken when q — q asymmetry is broken and 
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FIG. 4. (color online) v 2 of p and p at 7.7, 11.5 and 39 GeV 
for 0-80% calculated within UrQMD 



the width of distribution is sufficiently narrow. When 
a 2 = 0.25, the proton V2/n q is below the meson one by 
nearly 50% while the antiproton one is above it. K + 
and K~ arc similarly very different, while the difference 
between A and A is smaller, see Fig. 3, because the dif- 
ference comes only from one over three quarks instead 
of one over two. All of these can be readily explained 
by the greater admixture of high-flow antiquarks and 
strange quarks in anti-baryons and strange hadrons. In 
the K~(su), for instance, both quarks tend to come from 
smaller \y\ where the flow is larger, while in K + only 
one does. Hence, K + (u~s) has a higher V2(pr)- Up to 
an asymmetry of a y = 0.5 the breaking of the scaling is 
still sizable but start to emerge when the transverse mo- 
mentum per constituent quark is > 1 GeV. With a y > 1 
the usual scaling arc essentially restored. We notice that 
in a realistic description also thermal smearing should be 
included destroying the "locality" in momentum rapid- 
ity of the parton distribution and mixing quarks from 
flowing central region with the more baryon-rich frag- 
mentation region. This effect with Ay ~ 0.5, produces 
a similar effect to the one we have investigated for co- 
alescence with finite width wave function, therefore the 
lowering of baryon V2 would persist even if Aligner — > 
as long as p/p breaks boost invariance. 

Our other main result, showing in Fig. [31 is that V2 
is higher for anti-baryons than for baryons. The differ- 
ence is noticeable at all px if the width in rapidity of the 
antiquark and strange quark distribution is 0.25, but be- 
comes noticeable at higher px ~ 2 GeV only, if the width 
increases to ~ 0.5. A further increase of the width makes 
the limiting px of the difference even higher. 

This result is stable against changes in the details of 
parton distributions as long as Eq. [3] and the decreasing 
q/q as a function of rapidity are valid assumptions, since 
baryons admit an admixture of quarks from the high 
rapidity peripheral regions, where V2 is smaller, while 
anti-baryons are dominated by central rapidity region at 
larger V2 ■ Any coalescence incorporating finite width 
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wave function will result in a higher v% for anti-baryons 
w.r.t. baryons. However, this directly contradicts ex- 
perimental data [28|, suggesting either that coalescence 
breaks down at lower energies, or that quarks and anti- 
quarks do not have the same flow i>2(j/,Pt)- 



IV. DISCUSSION 

At this point one may ask if at lower energy what we 
have created is just a pure hadron gas. The systematics of 
u 2 (pt): overlapping for all accessible momenta between 
yfs of 7.7 GeV and LHC energies (2^1, makes different 
regimes of v-i origin immediately suspect. Nevertheless, 
to investigate also this possibility we have calculated the 
V2 of protons and antiprotons in purely hadronic molec- 
ular dynamics, implemented via the UrQMD model 371 . 

It should be noted that the version we used, v2.3 38j . 
suffers, as in earlier versions, from lack of detailed bal- 
ance for mult i- particle reactions: annihilation processes 
involving multi-particle final states, such as pp —> iriririr 
are possible. The corresponding creation processes, how- 
ever, are absent. Quantitatively, this is a minor correc- 
tion, since its equilibration time for 2 -<-> n processes 
T n^2 ~ e" _2 r 2fl 2 1 and, as uRQMD shows, even 2 <-> 2 
processes do not equilibrate in realistic heavy ion collision 
expansion profiles. However, in high chemical potential 
systems antiprotons themselves are also minor correction, 
being suppressed by factors of ~ exp(— (m + (j,b)/T), so 
the relative importance of multi-particle processes could 
become enhanced. 

The uRQMD result is shown in Fig. \M for the upper 
and lower physically relevant energies [28| , it has a simi- 
lar behavior to purely partonic dynamics+coalescence. 

The reason this happens is however somewhat dif- 
ferent: Antiprotons are produced initially, and get ab- 
sorbed by annihilation in a proton-rich medium. Neglect- 
ing regeneration, it is not surprising that in a hadronic 
medium protons and antiprotons propagate differently: 
Their thermally averaged interaction cross-sections are 
considerably different, since protons typically interact via 
resonance formation while antiprotons, in a high chemical 
potential medium, can interact both via resonance forma- 
tion and collisional processes and annihilation, with the 
latter dominating if the Knudscn number is large and V2 
is created by absorption rather than flow. This, how- 
ever, already makes it likely that antiprotons be more 
anisotropic, both due to flow and absorption, simply be- 
cause they interact more with the medium. 

Once again, naive extrapolation seems to disagree with 
experimental data [28], either purely hadronic dynam- 
ics nor the coalescence hadronization of a partonic fluid. 
Ref.(30j proposed to explain this by assuming V2 of par- 
ticles transported from high rapidity is higher than the 
i>2 of particles created at mid-rapidity. It is however dif- 
ficult to see how such an ansatz is compatible with the 
scaling of v- 2 with rapidity [H, HJ, since "transported" 
degrees of freedom have to come, and spend more time 



in, less flowing regions. Our work confirms, in fact, that 
both partonic and hadronic dynamics, considering longi- 
tudinal diffusion only, have the opposite dependence on 
rapidity than that claimed in [30j . 

As shown in (39j . transverse ideal hydrodynamics, in- 
corporating instant chemical and thermal equilibrium for 
all hadronic species, might provide a way out, via chem- 
ical composition anisotropy in the transverse plane. If 
flow is created close to the ideal hydrodynamic limit, 
outer regions in transverse space have significantly more 
flow than inner-lying regions. The p/p ~ e ~ 2 ^B/T ra ^j 0; 
however, decreases with transverse radius since most p 
are formed in the higher density inner regions: in the 
outer regions \ib is approximately the same but T goes 
down, depleting the p/p ratio. 

As shown in [HI, this effect can drive i>2 of protons 
above i>2 of antiprotons. However, coalescence in this 
case is also not ruled out, since by exactly the same rea- 
soning, given a partonic system with q — q transverse 
anisotropy, one can drive the mean flow of quarks above 
the mean flow of anti-quarks, effective relaxing the as- 
sumption of Eq. 2J V2 q {y) = V2q(y): locally in transverse 
space flow is of course still the same for all species, but 
the inhomogeneity leads to different transverse-averaged 
flows. Fig. [5] shows what happens if the effective V2 q 
is just 20% of those of the antiquarks, a comparable 
effect to that used in [39| (the exact magnitude of the 
q — q transverse anisotropy should be approximately in- 
sensitive to whether uRQMD or "partonic" initial con- 
ditions are used). Partonic chemical non-equilibrium at 
freeze-out might accentuate this imbalance [401, since at 
fis 3> T light quarks equilibrate radiatively while light 
antiquarks equilibrate collisionally. 

We conclude, therefore, that 2+1 partonic hy- 
dro+coalescence might also be consistent with lower en- 
ergies RHIC data if the transverse radius dependence of 
elliptic flow is taken into account. Since, as discussed 
previously, strange quarks are chemically similar to light 
antiquarks in that they are created thermally, we expect, 
if this scenario is correct, that 



Vo > Vn 



Vo > Vn 



Alternatively, differential quark-antiquark flow can be 
directly realized in a system with either partonic or 
hadronic degrees of freedom and strong vector mean 
fields, since the repulsive channel at finite net quark num- 
ber will generally be greater for q, s than q, s (NB,both 
light and strange quarks). Vector mean fields, unlike 
scalar ones, can produce a difference between quarks and 
antiquarks because they admit both attractive and re- 
pulsive channels, and are sensitive to a conserved charge 
density such as the net baryon or quark density. 

In a medium with positive chemical potential, vector 
mean fields should cause relatively more attraction be- 
tween anti-particles and the higher baryonic density re- 
gions, and therefore a smaller V2(pt) for antiparticlcs 
with respect to the particles. As shown in [41| within 
Boltzmann transport theory an attractive mean field 
cause a significant reduction of the V2{pt)- This can 



6 



ultimately be related, in a hydrodynamical picture, to 
the fact that an attractive mean field decreases the effec- 
tive pressure of the system at equilibrium. While elliptic 
flow is thought to be driven by pressure gradients, it re- 
mains to be clarified up to what level the transport and 
hydrodynamic picture are equivalent. 

At hadronic level the vector mean fields are, ultimately 
how the AMPT calculation in [42j has been able to re- 
produce both the baryon and anti-baryon V2 data. 

A model with quark mean fields has also been recently 
conjectured and studied within NJL [4l[ and PNJL mod- 
els [43| . As piUI^ have shown quark degrees of freedom 
cou pled to mean fields can produce sizable elliptic flows 
in a vector channel, can again vary con- 
siderably between quarks and antiquarks. 

If, in the partonic system, the effect of mean fields is 
comparable to isotropic pressure, however, the very defi- 
nition of thermalization and r\j s needs to be revised, since 
defining a co-moving frame within a strong mean field is 
impossible: Particles and antiparticles will always move 
in different directions within a volume element, exactly 
the effect required to explain energy scan data. This also 
means that if V2 is driven by mean fields, one would ex- 
pect the difference between baryon and anti-baryon v-i 
to be constant with strangeness, since QCD dynamics is 
flavor-blind 

v\ - v p 2 ~ v£ - ~ uf - uf ~ v$ - v$ 
This latter scenario, however, leads to some unsettling 
questions for the interpretation of w 2 at higher ener- 
gies: Since at p/p ~ 1 such effects become invisible (C- 
symmetry of the medium is restored), there is no indica- 
tion that they actually turn off at RHIC/LHC, and the 



smoothness of the decrease of the 



with 



,/s [28| 

suggests they do not. As mean fields are "large" devia- 
tions from the hydrodynamic limit, discovering that they 
have a contribution to creating (rather than suppressing) 
v n could significantly modify viscosity estimates such as 
those at higher energies @. 

To constrain the models further using existing data, 
let us consider that the total (n, K,p) v^Ipt) in the same 
energy scan can be experimentally shown [29( to be ap- 
proximately independent of the initial density (which de- 



pends, approximately, as 



Npart m - 



1 



assum- 



ing Bjorkcn expansion [33( with a constant starting time); 
The integrated vi varies because (px) increases system- 
atically with -jt^, but V2(pt) overlaps within error bar 

& 

The analysis in Fig. |4] also fails to reproduce such a 
scaling specifically for protons, in line with the interpre- 
tation of absorption as the dominant origin of vi ; In this 
case, the V2 scaling with pr should be governed by the 
collision integral formula of [48j , where transverse density 
enters both into (pr) and V2{pt) via 



V 2 (Pt) ~ ((TVijipT)) 



e dN 
S dy 



(15) 



negative velocity for that particular transverse momen- 
tum. Neglecting averaging over rapidity, and assuming 
daij ~ aijdp^j/Q 4 where Q is the momentum transfer 
and ptj is the momentum of the second scatterer, we get 



/ dN 
v 2 {pt) ~ e-W 
dy 



Pt \ G(p T ) 
(pt)J S 



dc> 



f 



PTJ 



(16) 
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(Pt) J (p T + 



Ptj 



2ptPtj cos(0))^ 



where we used the universal scaling of the transverse 
momentum distribution with pxj (pt) noted in (49| and 
/(...) is a distribution function, which can be well ap- 
proximated by a Tsallis spectrum [50j . 

The G(pt) term generally breaks the scaling, increas- 
ing V2(j>t) by virtue of an increase in (pr) with y/s, N par t 
and decreasing v 2 (j>t) at the high pr ^> (pt) tail, 
in case of high baryochemical potential limit at low 
a/s, the averaging of (pr) over p,ir,K,... and ay over 
p — p,p — ir, tt — 7r, ... further breaks the scaling. This 
systematic shift is also visible in Fig. 11 of [29| or in [5l[ 
within the uRQMD and AMPT models, as v 2 (j>t) in- 
creases systematically with -y/s. The same reasoning ex- 
plains the systematic rise specific to p, p observed in Fig. 

H 

However, as explicitly discussed in [52, [H[ but clear 
from 0, [54[, an ideal fluid stage creating v 2 via az- 
imuthally inhomogeneus pressure gradients will not ame- 
liorate this problem, since a Cooper-Frye freeze-out [55| 
will result in a i>2 which to leading order depends on the 
2nd Fourier components of the freeze-out hyper-surface 
£2 T and transverse flow e^ T via 



v 2 (pt) ~ ytanh 



l e 2 Pt 



T 



(17) 

is done over events at a given ^/s 



Here (o~Vij(pT)) is the normalized cross-section times 



and the averaging (. 
and centrality class. 

Since both T and t v 2 depend on the lifetime of the 
system (albeit approximately saturate after a finite time 
[56(), V2(pt) should also systematically rise with y/s. The 
systematic deviation should increase with pt , since both 
e^ T '" T come with factors px- Such a behavior is indeed 
evident in Q , and is qualitatively the same as the trans- 
port simulation [5lj |. 

In the recent years it has been point out that an agree- 
ment with data for several observables even at the highest 
RHIC energy can be better explained if one introduce a 
core-corona model [13] ■ The fact that AMPT, which in- 
corporates both hadronic and partonic components in the 
relevant energy density regimes, has the same qualitative 
dependence as Q (as seen in [29| Fig. 11) indicates that 
a core-corona model might find it challenging to repro- 
duce V2(pt) for all energies, since the greater preponder- 
ance of the corona at lower y^i should also change the 
Pt dependence of V2 («2 in the corona is either zero or 
hadronic). However, since a core-corona model with an 
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P T [GeV] 

FIG. 5. (color online) Plot of baryon vs anti-baryons v% with 
same parameters as in Fig[3] (upper panel) but with a modified 
width for quark, see text for details 

ideal core (77/s = 0) does a good job of reproducing inte- 
grated v 2 /e at several top energies and centralities (57j . 
and n/s = is expected to scale the best [20[, perhaps 
the scaling of v 2 will to a certain extent survive when 
binned in px- This can be expected also on the base of 
dynamical transport calculations that include the strong 
Tj/s increase at hadronic energy densities (58j . Regard- 
ing the baryon/ anti-baryon pattern of V2(j>t), we have 
shown in the present paper that both a partonic and 
a hadronic picture show the same baryon/ anti-baryon 
pattern. Therefore we can say that also employing a 
core-corona picture the issue raised might persist unless 
strong transverse gradients within the core, a la (39j , are 
maintained. It is therefore difficult to see whether cither 
a hybrid analysis of a uRQMD-only analysis can accom- 



modate the difference between p and p as well as the 
near-overlap of V2(pt) seen in [29j. This near-overlap has 
not, however, as yet been established for separate particle 
species, so perhaps baryonic v 2 is quantitatively different 
from 7r-dominatcd total v 2 . In this case, hydrodynamic 
evolution with a long hadronic afterburner [39( might be 
the simplest explanation of elliptic flow systcmatics. 

In conclusion, we used quark coalescence and a sim- 
ple parametrization of a quark fireball with hydrody- 
namic flow to investigate the dependence of coalescence- 
type scaling at various energies. We concluded that it is 
wrong to expect coalescence to give the naive quark num- 
ber scaling of v 2 at all energies and system sizes even if 
partonic hydrodynamics and coalescence actually occur. 
We have also found that both coalescence and hadronic 
dynamics predict the opposite difference of baryon and 
anti-baryon v 2 w.r.t. experimental data. The scaling of 
baryonic v 2 , therefore, points to the existence of separate 
quark and antiquark flow at lower energies. A theoretical 
explanation of this will be essential to link flow dynamics 
in heavy ion collisions to the deconfinement phase tran- 
sition and this is presently missing. 
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